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ABSTRACT 


We have obtained ALMA Band 7 observations of the FU Ori outburst system 
at 0.6" xO.5" resolution to measure the link between the inner disk instability 
and the outer disk through sub-mm continuum and molecular line observations. 
Our observations detect continuum emission which can be well modeled by two 
unresolved sources located at the position of each binary component. The in¬ 
terferometric observations recover the entire flux reported in previous single-dish 
studies, ruling out the presence of a large envelope. Assuming that the dust 
is optically thin, we derive disk dust masses of 2 x 10 _4 M o and 8 x 10 -5 M 0 , 
for the north and south components respectively. We place limits on the disks’ 
radii of r <45 AU. We report the detection of molecular emission from 12 CO(3- 
2), HCO + (4-3) and from HCN(4-3). The 12 CO appears widespread across the 
two binary components, and is slightly more extended than the continuum emis¬ 
sion. The denser gas tracer HCO + peaks close to the position of the southern 
binary component, while HCN appears peaked at the position of the northern 
component. This suggests that the southern binary component is embedded 
in denser molecular material, consistent with previous studies that indicate a 
heavily reddened object. At this angular resolution any interaction between the 
two unresolved disk components cannot be disentangled. Higher resolution im¬ 
ages are vital to understanding the process of star formation via rapid accretion 
FU Ori-type episodes. 


Subject headings: FU Ori - Star Formation - circumstellar matter 
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1. Introduction 


It is well understood that stars acquire most of their material through an accretion 
disk that is fed, at least at early times, by a surrounding envelope. Once the bulk of the 
envelope is dissipated, the accretion rates tend to approach 10 7 or 10 8 M 0 yr h A 
small percentage of young stars - FU Orionis or FU Ori objects - exhibit rapid optical 


brightening of several magnitudes over timescales as short as one year (Herbig 1966). 
During the outburst phase, the accretion rates increase by factors of 10 3 or more relative to 
the quiescent periods, reaching rates as high as 10 -4 M Q yr -1 . This high state can last for 
~100 years and deposit a significant amount of disk material (~0.01 M 0 ) onto the young 


star in a short period (Hartmann & Kenyon 1996, and references therein). 


The cause of these outbursts is not well understood. Different mechanisms have been 
proposed as possible origin for the outburst (Audard et al. 2014): 1) thermal instability (Bell 


& Lin 1994), 2) coupling of gravitational and magnetorotational instabilities (Armitage et 


at] 2001; Zhu et al. 2009), 3) disk fragmentation (Vorobyov & Basu 2006; Zhu et al. 2012), 


and 4) tidal interactions between the disk and either planets (Lodato & Clarke 2004) or 


stellar companions (Bonnell & Bastien 1992). 


Connecting infrared studies which probe inside radii of 1-10 AU to the main mass 
reservoirs in the outer disk and/or envelope is particularly important if a full understanding 
of the phenomenon is to be achieved. For instance, it remains unclear whether the reservoir 
of material feeding the enhanced accretion rate is distributed in a circumstellar disk, in a 
larger envelope, or both. The only way to accurately recover the surface density distribution 
of the gas and dust around FU Ori objects is via resolved imaging of the disk at, preferably 
multiple, (sub)millimeter wavelengths and/or transitions where the emission is most likely 
optically thin. 


If the disk is small, and if there is no reservoir of material from which to replenish it, if 
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it continues to accrete at 10 _4 M o yr _1 , the disk may be quickly dying. If the disk is indeed 
more massive and/or larger than previous submillimeter observations indicate, it is then a 
counterpart to the more massive disks seen in nearer star forming regions and we are simply 
seeing one outburst of potentially many. Clearly, resolving the source and determining the 
mass is key. 


FU Ori 


Despite being the prototype of the FU Ori class and one of the nearest such objects, at 


a distance of 450 pc (Murdin & Penston 1977), little is known about the global properties 


of the accretion disk surrounding FU Ori itself. The estimated age of the system is ~2 Myr 
( Beck fe Aspin||2012 ). Early interferometric studies at 1.3 mm by McMuldroch (1995) found 
an unresolved (at ~3.5"resolution, or ~1500 AU ) disk with 13 CO(2-l) emission spanning 
8.1-13.3 km s -1 with a total integrated intensity of 1.4 Jy. Studies with the JCMT by 


Sandell & Weintraub (2001) show an unresolved continuum source (<5"in diameter) and a 
flux of 66 and 400 mJy at 850 and 450 pm, respectively. These results suggest there is no 
envelope bigger than 1500 AU around the source. 

Significant progress has been made towards understanding the inner disk of FU Ori 


objects in general and FU Ori specifically (see papers by e.g. Zhu et al. 2009; Eisner 


& Hillenbrand 2011). The region of rapid accretion for FU Ori itself has been shown 
to be occurring over a radius of ~1AU, or 0.005", in diameter. Given the size of the 
region involved, thermal instability alone cannot explain the outburst and appeals to other 
outbursts mechanisms must be made. 


FU Ori was discovered to be a binary system by (Wang et al. 2004), with the newly 


discovered binary component (FU Ori South, or FU Ori S hereafter) located 0.5"to the 
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South of the brighter principal component, FU Ori. Recent studies indicate FU Ori S 


is actually the more massive star in the system (~1.2M Q ; Beck & Aspin 2012), heavily 
embedded in obscuring material (Pueyo et al. 2012). A number of FU Ori type objects have 
now been identified as having binary components (e.g., L1551 IRS5, RNO 1B/C, AR 6A/B 


Pueyo et al. 2012, and references therein), triggering discussion on whether the rapidly 
accreting FU Ori phenomenon could be caused by close-companion interaction. However, 
at least in the case of FU Ori itself, this hypothesis remains inconclusive (e.g. see the 


discussion in Beck & Aspin 2012). 


Sandell & Weintraub (2001) suggest an upper limit for the disk mass at 0.02 M e . There 


are a variety of ways to constraint the dust, gas and total mass of the disk, all of which 
involve obtaining resolved imaging to determine the optically thick contribution and model 


fits (Andrews & Williams 2007; Eisner et al. 2008). Below we present new observations of 


FU Ori with sufficient resolution to constrain the masses of disks around FU Ori N and S. 


2. Observations and data reduction 

FU Ori Band 7 observations were acquired on December 2 nd 2012 using 24 antennas. 
The dataset provides baselines from 15 to 381 meters. During the observations the 
precipitable water vapor in the atmosphere was stable between 1.3 and 1.7 mm with clear 
sky conditions, resulting in median system temperature of 190 K. The ALMA correlator 
was configured in the Frequency Division Mode (FDM) to provide 468.75 MHz bandwidth 
in each of four different spectral windows at 122 kHz (~0.1 km s -1 ) resolution per channel. 
Three spectral windows were positioned to target the 12 CO J = 3 — 2, HCN J = 4 — 3 
and HCO + J = 4 — 3 transitions at 345.79599, 354.50547 and 356.73424 GHz, respectively. 
The fourth spectral window was positioned in a region devoid of line emission for detecting 
the continuum dust emission (centered at 345.78359 GHz). Callisto was observed as flux 
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calibrator, while the quasars J0522-364 and J0532+075 were observed for bandpass and 
phase calibration respectively (J0532+075 is located 3.5 degrees away from the science 
target). Observations of the phase calibrator were alternated with the science target every 
5 minutes to calibrate the time-dependent variations of the complex gains. The total time 
spent on-source was 16.9 minutes. A secondary phase calibrator, J053056+13322, was also 
observed regularly (located 5.7 degrees away). 


All the data were calibrated using the Common Astronomy Software Applications 
package (CASAQ |McMullin et al. 2007) in a standard fashion, which included offline 
Water Vapor Radiometer (WVR) calibration, system temperature correction, as well as 
bandpass, phase and amplitude calibrations. After flagging for problematic antennas, a 
total of 21 antennas were used for imaging. The fluxes derived for J0522-364, J0532+075 
and J053056+133220 were 7.5, 1.1 and 0.5 Jy respectively. During Cycle 0, the absolute 
flux calibration of the data is estimated to be accurate to within 10 to 15%. 


Imaging of the continuum and molecular emission lines was performed using the 
CLEAN task in CASA. Using natural weighting resulted in a synthesized beam size of 
0.63" x0.51" at PA= 3.0 degrees. Continuum subtraction in the visibility domain was 
performed prior to imaging of each molecular line. After CLEANing the images, an RMS 
noise level of 25.0 mJybeam -1 per 0.22 km s -1 channel was reached (in the 12 CO J = 3 — 2 
line). CLEANing of the dust continuum was performed after combining the line-free 
channels from all four spectral windows (adding to a total bandwidth of 1.759 GHz centered 
at 351.33 GHz), for which an RMS of 0.64 mJybeam -1 was obtained. After a single 
iteration of phase-only self-calibration the RMS of the continuum image was reduced to 
0.31 mJybeam -1 . 

J http://casa.nrao.edu/ 
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3. Results 


3.1. Dust Continuum 


Figure [l] shows the ALMA Band 7 continu um image of FU Ori after self-calibration. 
The continuum emission at 0.854 mm (351.33 GHz) is detected at a signal-to-noise ratio 
(SNR) of ~150. The emission is resolved, with the peak roughly coinciding with the 
position of the northern binary component (FU Ori). We used the CASA task IMFIT to fit 
the observed emission using a single component 2-D Gaussian, but inspection of the residual 
image shows that a single-component fit fails to match data. Using two unresolved sources 
provides a much better fit to the observed emission (with residuals down to the image 
noise levels). The total fluxes of each component are 50.1±0.3 mJy and 21.2T0.4 mJy 
respectively. 


The parameters of the Gaussian components determined using IMFIT are detailed in 
Table [lj The positions of each Gaussian component are shown in Figure [!} The position 
of the brighter component coincides within 0.1” of the optical position of FU Ori (as listed 
in 2MASS), whilst the fainter disk is located 0.5” to the South-East (0.48” to the South 
and 0.14 to the East). This is similar to the separation between the FU Ori and FU Ori S, 


as reported in previous near-IR observations (e.g. Wang et al. 2004; Monnier et al. 2009 


Beck & Aspin 2012; Pueyo et al. 2012, and references therein). 


3.2. Line emission 

We detect emission from 12 CO J = 3 — 2, HCO + J = 4 — 3 and HCN J = 4 — 3 
molecules at SNR levels (peak/rms) of 47, 12 and 5 respectively. Figure [I] and Figure [2] 
show the integrated emission maps for each molecule. The integrated line emissions are 
31.61±0.12, 2.14±0.07 and 0.37±0.07 Jy km s -1 respectively. Figure [3] shows the 12 CO 
















channel maps towards FU Ori. The integrated intensity maps of the three molecules are 
strikingly different. 12 CO is more extended compared to the continuum (extending up 1” 
from FU Ori, and covering the two binary components), and peaks east-ward from FU Ori 
On the other hand, HCO + peaks close to the position of the southern binary component 
(FU Ori S), whilst HCN J = 4 — 3 appears to peak closer to FU Ori. 


Figure [ 3 ] shows the 12 CO channel maps towards FU Ori. There is significant cloud 
contamination in the central channels, making the interpretation of the 12 CO spectrum 
difficult. The velocity of the channels with no 12 CO emission (11.8 km s -1 ) is very similar 
to the velocity channels where strong 12 CO(1-0) emission was previously detected with the 
NRAO 11m telescope (Kutner et al. 1982). This suggests that unresolved emission and/or 
absorption from the larger Orion molecular clouds contaminates the emission in the central 
channels of the ALMA spectrum. Nevertheless, a velocity gradient can be perceived in the 
channel maps and in the first moment map of the 12 CO (Figure [4]). Figure [ 4 ] shows the 
position-velocity diagram taken in the East-West direction (where the velocity gradient 
appears to be larger), which suggests there is gas in keplerian rotation around the system. 


4. Discussion 


We detect continuum emission towards both binary components, suggesting they both 
harbour circumstellar disks. The disks are unresolved at the M).6"x0.5" resolution of our 
observations (or 225 AU at FU Ori’s 450 pc distance), placing upper limits on the disk 
radial sizes. To be consistent with the Gaussian fit, the disks must have radial sizes smaller 
than 0.1” (45 AU). At the projected separation of the binary system (225 AU), the disks 
radii are predicted to be truncated at 0.2 to 0.5 times the binary separation (Artymowicz 


& Lubow 1996). The deconvolved disk sizes of 45 AU coincide with 0.2 times the binary 


separation. 
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Our ALMA data rule out the presence of a large envelope around the two stars, since 
our observations would have detected any extended emission smaller than 7” (3000 AU) 
in diameter. Our observations also place limits on the presence of a possible circumbinary 
ring orbiting the system. Assuming a 0.5” wide ring, starting at 1.5” (given the minimal 
truncation radius for the projected binary orbit; e.g. Artymowicz & Lubow 1996), the ring 
would need to have a dust mass smaller than 1.9 x 1O _4 M 0 to remain undetected by our 
observations (at 3-<r, assuming a dust temperature of 20 K). This is similar to recent ALMA 


results from Akeson & Jensen (2014) in which no circumbinary emission was detected in 


any of the 17 binary stars of their sample. 

Assuming a dust temperature of 50 K for the individual disks around FU Ori and 


FU Ori S, a gas-to-dust ratio of 100, and a 345 GHz dust opacity of 1.5 cm 2 g 


2 „-l 


we can 


compute a coarse estimate of the disks’ total masses (e.g. Beckwith et al. 1990). Under 
these assumptions, the ALMA 854 pm fluxes correspond to disk dust masses of 2 x 10 -4 M 0 
and 8 x 10 -5 M 0 , for the north and south components respectively. The total mass of the 


combined disks is comparable to the mass estimated from unresolved observations (Sandell 


& Weintraub 2001). These mass estimates are, however, subject to the assumption that 


the dust emission is optically thin. For the disks to be optically thick, the average disk 
dust temperature would have to be higher than 80 K. Assuming a R~ 0 5 dependence for 
the disk temperature, an average T > 80 K is reached if the temperature at 45 AU is 
greater than 60 K. This is high for a typical T Tauri disk, but it corresponds roughly to the 
radiative equilibrium temperature for a source with a luminosity larger than about 4 L 0 . 
This is perfectly feasible for FU Ori, where the outburst luminosity much larger than the 
photospheric luminosity. We conclude that the disks could be marginally optically thick, 
although higher resolution imaging would be required to accurately constrain this. 


According to Beck & Aspin (2012) the masses of FU Ori and FU Ori S are 0.3 M 0 , 
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and 1.2 M 0 respectively. Assuming an effective temperature of 4000-6500 K and an optical 


extinction of A\/=8-12, Pueyo et al. (2012) estimate a minimum mass for FU Ori S of 


0.5 M 0 . Both studies agree that FU Ori S is the most massive star in the system. In this 
scenario FU Ori S would be the primary star in the system. Unfortunately, our line data 
lacks the angular resolution to resolve the velocity gradients across each disk, which could 
be used to confirm these stellar masses. 

Theoretical models predict that the circumprimary disks should be the most massive 


disks in the system (e.g. Bate 2000). However, recent submillimeter surveys of binary 


systems show no clear relationship between disk masses and stellar mass ratios (Harris 


et al. 2012; Akeson & Jensen 2014). As noted by these authors, the predictions from 


these theoretical models only consider the initial conditions of the binary disks and do not 
take into account the dynamical evolution of the disks (e.g. viscous accretion and other 
dissipative processes). Assuming that both disks have the same average dust temperature, 
we find that in the FU Ori system the brightest and most massive disk orbits the less 
massive star. We note that this scenario may be different if the southern disk is cooler than 
the northern disk. For instance, if the southern disk has an average dust temperature of 
20 K, its mass would be higher than the 2 x 10 -4 M 0 dust mass derived for the northern 
disk. 

The denser gas tracer HCO + J = 4 — 3 peaks in the position of the most massive star 
FU Ori S, and is not detected in the direction of FU Ori. Since HCO + is a higher density 
gas tracer, this suggests that FU Ori total disk mass is probably lower than FU Ori S’s 
disk. On the other hand, HCN is detected only toward FU Ori. HCN forms at large heights 
above the midplane, where the temperature are higher and the densities are lower compared 


to the midplane (e.g. Walsh et al. 2010). Therefore, the HCN detection around FU Ori 


and not FU Ori S’s could support the idea that the disk around FU Ori has a higher dust 
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temperature compared to the disk around FU Ori S. As noted above, an average dust 
temperature below 20 K for FU Ori S would make it the most massive of the two disks, 
which would be in agreement with the interpretation of the HCO + detection. 


Based on the stellar photospheric absorption features Beck & Aspin (2012) estimate 
the accretion rate of FU Ori S in (2 — 3) x 10 _8 M o yr _1 . For FU Ori the accretion rate 
derived by Zhu et al. (2007) is 2 x 10 _4 M Q yr _1 . However this is highly model dependent 
is likely to correspond to the maximum accretion rate during the peak periods of episodic 
accretion. It is unlikely that FU Ori’s high accretion rate can continue for long. If FU Ori 
continues to accrete at 10 -4 M 0 yr _1 , its submillimeter flux would decrease by a factor of 
ten within a decade. To sustain the disk dust mass we infer from our ALMA observations, 
the stellar accretion rate must be highly variable. 


5. Conclusions 

We obtained ALMA Band 7 observations of the FU Ori binary system at 0.6" x 0.5" 
resolution in Cycle 0. Our observations detect continuum emission from two unresolved 
sources at the position of each binary component. This indicates that both binary 
components have disks. Under simple assumptions, these observations suggest that the 
most massive disk orbits the less massive star (FU Ori). 

The interferometric observations recover the entire flux reported in previous single-dish 
studies, ruling out the presence of a large envelope for replenishing the disks. The dust 
disks were unresolved in our observations, implying that they are significantly smaller than 
non FU Ori objects such as e.g. AB Aur, HL Tau and HD 163296. 

We report the detection of molecular emission from 12 CO(3-2), HCO + (4-3) and from 
HCN(4-3). Whilst the 12 CO appears widespread across the two binary components , HCO + 
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peaks close to the position of the southern binary component. This suggests that the 
southern binary component is embedded in dense molecular material, which is consistent 
with infrared studies that indicate that this is a heavily reddened object. Is this clump 
the reservoir of material feeding the outburst of the northern binary component? At this 
angular resolution the interaction between the two unresolved disk cannot be disentangled 

Observations at higher spatial resolution are mandatory to investigate whether 
interactions between the disks or other mechanisms (e.g. disk fragmentation) are causing 
the outburst. 
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Table 1: Integrated Fluxes and Gaussian Fit Parameters 


Component 

Gaussian 1 

Gaussian 2 

Position (Ra, Dec) 

(05:45:22.36,+09.04.12.24) 

(05:45:22.37,+09.04.11.75) 

Major axis (mas) 

662.3 ± 3.4 

636.3 ± 9.0 

Minor axis (mas) 

536.8+4.3 

517.9 ± 8.9 

Position angle (deg) 

0.8 ± 0.7 

9.9 ± 2.3 

Peak Intensity (mJybeam -1 ) 

45.2+0.3 

20.9 ± 0.3 

Integrated Flux (mJy) 

50.1+ 0.3 

21.2 ± 0.4 



Fig. 1.— Left-panel: ALMA Band 7 continuum image of the FU Ori system with contours 
overlayed (after self-calibration). Contour levels are 10, 20, 40 and 80x0.31 mJybeam -1 
(the rms noise). The stellar position from 2MASS is shown with a black star symbol. 
Right-panel: 12 CO(3-2) integrated intensity (moment 0) image with Band 7 continuum con¬ 
tours overlayed (Contour levels are identical in both figures). The peak 12 CO(3-2) flux is 
5.8 Jy beam -1 km s -1 . 
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RA (J2000) RA 02000) 


Fig. 2.— Left-panel: HCO + (4-3) integrated intensity (moment 0) image with Band 7 contin¬ 
uum contours overlayed. The peak flux is 0.88 Jy beam -1 km s -1 . Right-panel: HCN(4-3) 
moment 0 image with Band 7 continuum contours overlayed. Contour levels are identical to 
Figure [I] The peak flux is 0.37 Jy beam -1 km s -1 . 
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Fig. 3.— CO channel maps towards FU Ori. The velocity of the channels is shown in the 
Local Standard of Rest (LSR) frame, centered at the rest frequency of 12 CO J = 3 — 2. 
The data has been binned to a velocity resolution of 0.21 km s -1 . The band 7 continuum 
contours are overlayed (with contour levels are identical to Figure [I]). Cloud contamination 
is present in the channels between 11.63 and 12.69 km s -1 . All maps share the same linear 
color scale, ranging from -0.2 to 1.5 Jybeam -1 . 
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Fig. 4.— Left-panel: CO moment 1 (intensity-weighted velocity field) of FU Ori. The 
position of the two fitted continuum point sources and the band 7 continuum contours 
are overlayed (with contour levels are identical to Figure [TJ. Right-panel: FU Ori position- 
velocity diagram for the observed 12 CO line emission, centered at the velocity of 11.8 km s -1 . 
The spectra ±1” on either side of the apparent major axis have been averaged to form this 
PV diagrams. Spectral/spatial resolution is shown in the lower left. For comparison, the 
solid lines show the maximum Keplerian velocity for gas at the tangential point (at three 
different inclination angles), as a function of distance from the star. The mass of the central 
object has been assumed to be 1.5 times the mass of the Sun (i.e. the mass of the combined 
binary system). Green, red and magenta correspond to inclination angles of 30, 45 and 
70 degrees respectively. 









